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ABSTRACT 

We report on the detailed analysis of the high-energy extended emission from the 
short Gamma-Ray Burst (GRB) 081024B, detected by the Fermi Gamma-ray Space 
Telescope. Historically, this represents the first clear detection of temporal extended 
emission from a short GRB. The light curve observed by the Fermi Gamma-ray Burst 
Monitor lasts approximately 0.8 seconds whereas the emission in the Fermi Large Area 
Telescope lasts for about 3 seconds. Evidence of longer lasting high-energy emission 
associated with long bursts has been already reported by previous experiments. Our 
observations, together with the earlier reported study of the bright short GRB 090510, 
indicate similarities in the high-energy emission of short and long GRBs and open the 
path to new interpretations. 

Subject headings: gamma rays: bursts 



Introduction 



Gamma-Ray Bursts (GRBs) are extremely energetic and brief explosions originating at cosmo- 
logical distances. Since their discovery in the soft gamma-ray regime roughly 40 years ago, they have 
been detected to emit in almost every wavelength, from mm to the GeV range. The properties of 
their gamma-ray emission (e.g., duration, spectral shape, variability) have been thoroughly studied 
in the past 20 years, with multiple spacecrafts, notably with the Burst And Tran sient Source Exper- 
iment {BATSE) onboard the Compton Gamma- Ray Observatory [CGRO) (see lZhang &: Meszaros 
20041 . fo r a review). One of the ir most enduring properties is their classification in two duration 
classes ( Kouveliotou et al. 1993 ). with distinct spectral characteristics: short ( ~ 2 s) hard GRBs 



and l ong soft ones. This bimodality has been confirmed with NA SA's Swift satellite (iGehrels et al 



20041 ) and recently with the Gamma-ray Burst Monitor ( GBM, iMeegan et al.l |2009| ) onboard the 
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Fermi satellite. The prevailing notion is that these two GRB classes originate from different pro- 
genitor systems. 

The first dete ction of an X-ray afterglow goes back to observation of GRB 970217 with the 
BeppoSax satellite (ICosta et al.lll997l ) , followed by the discovery of the optical tr ansient that allowed 
the first ever determination of the redshift for a GRB (jvan Paradijs et al.lll997l ). Multi- wavelength 
(radio to X-ray) follow-up observations of GRB afterglows and the extensive study of the GRB 
hosts and their environments, have now established that at least so me long GRBs are conn ected 
with the collapse of massive rapidly rotating stars into black holes (jWooslev &: BloomI |2006| . and 
references therein). The origin of short GRBs is less certain. In the last few years growing evidence 
supports the idea that short GRBs originate from neutron star-neutron star (NS-NS) mergers also 
ending in a few so l ar mass black-hole surrounded by a short-l i ved accretion disc (s ee 



19891 : 



RufFert fc Janka 



1999 



Rosswog et al.ll2000l : iNarayan et al 



2001 



Nakai 



2007- 



Rosswog et al 



Eichler et al. 

2OO3I ). It is. however, rather difficult observationally to distinguish mergers from collapsars. SN 



explosions have been spectroscopically associated in a scant four cases with GRBs, but none of 
these has been a short GRB. 

The detection of high-energy (^100 MeV) prompt and afterglow emission from GRBs with 
Fermi offers a unique probe to test the properties of the outflow, determine the mechanism 
of energy transfer, and to address fundamental p hysics issues like Lorentz invariance, a s has 



been done for other bursts such as GRB 080916C ( Abdo et al 



2009c 



Zhang fc Pe'eil l2009l ) and 



GRB 090510 (jAbdo et al.l l2009dl'l . During its first year of operati on Fermi detected high-energy 



photons from two short GRBs (jOmodei 



2008 



Qhno &: Pelassall2009l ). opening a new window in our 



understanding of the GRB phenomenon. 



High-energy emission above 1 GeV had been detected in the past with CG-RO/EGRET (jSommer et al 



1994 ) in association with several long BATSE GRBs. Little was known, however, about the ra- 
diation physics of short GRBs, due to poor statistics a t higher energies (IKaneko et al.ll2006l ). A 
notable exception is the extremely bright GRB 930131 ( Kouveliotou et al.lll994l ). which consisted 
of an initial very intense complex of pulses with total duration about 1.5 second, superimposed on 
a significant shallow tail lasting up to ~50 s. As the tail fluence was at least an order of magnitude 
lower than the one of the initial complex, it was the initial short pulse that dominated the event 
energetics and most likely contained most of the burst flux. Apart from its intensity, GRB 930131 
was also remarkable in its high-energy properties; EGRET detected 16 high-energy gamma rays, 
including the highest energy photon of 1.2 GeV 26.5 s after the initial spike. As the event was not 
classified as a short GRB, however, it was not clear until recently that short bursts could produce 
such high-energy emission, and whether the properties of the short GRB emission resembled those 
of the long events at GeV energies. 

In this paper, we report on the observation and the analysis of the GeV emission from the 
short GRB 08102 4B. The event lasts appr oximately 0.8 s below 5 MeV, entirely consistent with the 
short GRB class (jKouveliotou et al.lll993l ). The emission above 100 MeV lasts for about 3 seconds. 
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thus showing unambiguous evidence for a delayed high-energy component in a short GRB. 



2. Observations 



The Fermi Gamma-ray Space Telescope was launched on 2008 June 11, and during the first 
year of science operations detected high-energy emission from 9 GRb£|, doubling the number of 
bursts detected a bove 100 MeV. Two instruments operate onboard Fermi: the Gamma-ray Burst 
Monitor (GBM) (jMeegan et al.ll2009l) cover ing the energy range from 8 keV to 40 MeV, and the 
Large Area Telescope (LAT) (jAtwood et al.1 12009). from 20 MeV to more than 300 GeV. 



On 2008 October 24 at 21:22:41 (UT) the LAT detected an increased count rate associ- 
ated with the short burst GR B 081024B, which triggered the GBM (trigger number 246576161, 
(jConnaughton &: Briggsl l2008l )). The LAT ground an alysis for burst detection and localization 
follows the procedure described in lAbdo et alJ (j2009bl ) . Selecting "transient" events with energy 
above 100 MeV, final localization is found to be RA = 322°. 9 , Dec = 21°. 2 , with a statistical 
uncertainty of 0°.2 (68% containment radius). The map of the Test Statistics reaches its maximum 
value, TSmax = 48.9, at this position, corresponding to a 6.7 a d etection. We a lso estimated the 
significance of the high-energy emission using the Li &: Ma method (jLi &: Malll983l ) , which computes 
the probability of the temporal excess in the case of Poisson statistics. We used two different event 
selections, considering a fixed region of interest (ROI) of 15°, and an energy dependent ROI, taking 
into account the energy dependence of the LAT point spread function (PSF). The LAT recorded 
11 events with reconstructed energy above 100 MeV, within 15° from the position of the burst 
and within 3 s from the trigger time, when the expected number of counts from the background 
is 0.75, corresponding to an excess of 6.8 a. The number of expected counts has been computed 
considering the background 100 s before and 100 s after the burst, excluding the time window of 3 
s during the burst. A more careful estimation of the background, using Monte-Carlo simulations 
to estimate the charged particle background, and six months of data for estimating the gamma- 
ray backgro und was also performed, providing a significance of 6.7 a. Details of the method are 
described in lAbdo et al.l (|2009bl ). If we consider the energy dependent ROI, considering only the 
events that are within three times the 68% containment radius of the PSF (approximately corre- 
sponding to the 95% of the containment radius), the number of expected events from background 
decreases to 0.08, and the significance of the excess increases up to 8.5 a. This higher significance 
can be easily understood as the LAT PSF strongly depends on the energy (oc E~^'^), which makes 
an energy dependent ROI very effective in reducing the background contamination, without loosing 
gamma-ray events from the source. 

The LAT position was observed by the Swift X-Ray Telescope (XRT), startir ig 70.3 ks after 



the GRB trigger time and lasting 9.9 ks, but no X-ray counterpart was found ([Guidorzi et al 



' [http: //fermi .gsfc .iiasa.gov/ssc/resources/observations/grbs/grb_table7] 
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20081 ) ■ An additional 13.5 ks observation was conducted the next day (jGuidorzi &: Marguttil 120081 ) 
and confirm ed the lack of an X-ray afterglow candidate. Optical observations also produced no 
counterpart (|Cenko &: Kasliwa]||2008l ). The absence of an XRT detection so long after a short burst 
is not too surprising. A study of the Swift XRT catalog of GRbH shows that of the 33 short GRBs 
for which XRT observations were made, only 10 had detectable flux 75 ks after the trigger time. 



3. Light curves, durations and spectral lags 

The multi detector light curve is shown in Fig. [TJ The top panel shows the background 
subtracted light curve for the summed signal of the two brightest Nal detectors (6 and 9) between 
8 keV and 260 keV. The background subtracted light curve of the brightest BGO detector (1) is 
shown in the second panel (260 keV - 5 MeV). The third panel shows the LAT signal without 
any selection (i.e., all the events that passed the onboard gamma filter). The quality of these 
events is not good enough to use them in the spectral analysis, but the properties of the ensemble 
can be assessed quantitatively by means of a dedicated Monte-Carlo simulation and do convey 
physical information in the extremely interesting energy range in which the BGO and the LAT 
overlap. These three light curves are characterized by a narrow spike of about 0.1 s (interval "a"), 
followed by a longer pulse, of about 0.7 s (interval "b"). There is no evidence of emission after 
~0.8 s in the NaT and BGO detectors. The bottom panel shows the light curve of the "transient" 
selected events with well defined direction and energy (above 100 MeV). The arrival times and the 
reconstructed energies (right axis) of the selected events are also displayed. We estimate a total of 
0.4 background events during the time interval shown in Fig.[TJ An event with energy 3.1 it 0.2 GeV 
was detected after 0.55 s while a second event of 1.7 it 0.1 GeV was detected after 2.18 s. Tabled] 
contains the arrival times, the energies with the estimated error, and the arrival directions of these 
eleven selected events. The last two columns of the table are the estimated 68% containment radii 
calculated from the point-spread function (PSF) and the distance from the localization of the GRB, 
in PSF units. 

We studied the narrow spike visible in the full lightcurve (third panel) of interval "a". The 
probability to obtain the same number of counts from background fluctuations is discarded at the 
3.5 cr level. Furthermore, we performed dedicated Monte-Carlo simulations to estimate properly 
the energy of these LAT photons. These events do not belong to the "transient" class, which is the 
most generous event selection with minimal requirements on direction and energy reconstructions. 
Typically, they are discarded because they produce very few hits in the tracker (<20), with a very 
short track, and very low raw energy deposited in the calorimeter (<5 MeV). These topologies are 
typical of low-energy events, with energies between 10 MeV and 40 MeV. If we select these topologies 
in the data, the probability that this narrow pulse is the result of a background fluctuation decreases 
to 5.9 cr level. 



^http: //swift .gsfc .nasa.gov/docs/swif t/archive/grb_table 
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We conclude that the spike in interval "a" in the LAT data is significant, with energies below 
100 MeV (although the energy resolution for this class of events is relatively poor. ~50%, from 
Monte-Carlo simulations). 

A common method for estimating GRB durations is to compute the T90 which measures 
the du ration of the time inter val during which 90% of the total observed counts have been de- 
tected (jKouveliotou et al.lll993l ). Background fluctuations, especially in weak GRBs, strongly affect 
Tgo's as often the fluctuations are comparable to the 5% of the total GRB fluence. Values of T50, 
however, are generally more robust, and especially so in the case of weak events. For GRB 081024B, 
the T50 (T90) is 0.33 s (0.66 s) in the Nal detectors, 0.15 s (0.27 s) in the BGO, while it is signifi- 
cantly longer for the LAT, corresponding to 0.9 s (2.1 s) for the full statistic light curve and 1.5 s 
(2.6 s) selecting only the events above 100 MeV. 

Spectral lags are characteristic of long GRBs, which exhibit hard-to-so ft spectral evolution, 
while short GRBs do not exhibit such a property (jNorris &: Bonnelll l2006l ) . We searched for a 
possible spectral lag in GRB 081024B using the cross-correlation function (CCF). The GBM Time 
Tagged Event light curves of the 4 brightest Nal detectors (6, 7, 9 and 10) were summed with 
100 ms time resolution in 8 logarithmic energy bins from 8 to 1950 keV. Similarly, the two BGO 
light curves were added in 8 bins from 0.11 to 107.6 MeV. For the LAT, we us ed all photons 
above 100 MeV. The errors on the CCFs were estimated using Bartlett's formula (lBartlettlll978l ) 
and were propagated to the errors of the peak position. We computed the CCFs between the 
GBM/Nal and the GBM/BGO detectors as well as between the GBM and the LAT. We found no 
energy-dependent delay between any of the data types used. We also used two resolutions (50 and 
100 ms) and verified that our results did not change significantly. 



4. Spectral Analysis 

We performed a time-resolved spectral analysis in intervals "a" , "b" , and "c" of Fig. [H si- 
multaneously fitting the signal from the two Nal (selecting all the channels between 8 keV and 
860 keV), the BGO (from 200 keV to 36 MeV) and the LAT detectors (selecting "transient" events 
above 100 Mevjl. Table [2] shows the results for all time intervals, testing different fitting functions. 

In interval "a" the be st fit to the GBM data is obtained with a power-law with exponential 



cutoff (jKaneko et al.ll2006l . see COMPT model). The best fit parameters are summarized in Table[2] 
as fit "1" . The peak energy -Epeak lies in the BGO energy range, and, even though its value is 
only marginally constrained, is consistent with a very hard spectrum, with a roll-off at energies 
above a few MeV. The LAT upper limit on the photon flux in the 100 MeV-10 GeV energy range 



^Here we have used the post-launch Pass 6 v3 IRFs, which take into account, on average, the correction of the 
inefficiency due to pile-up of cosmic rays. The version of Science Tools used is vl5r9p2. We have used the package 
rmfit (version 3.1) to perform the GBM/LAT joint fit. 
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is 3.8 X 10~^ ph cm~^ or, in energy flux, 4.7 x 10~^ erg cm~^ s~"^, and is consistent with the 
extrapolated flux fr om the COMPT fu nction fitted to the GBM data. We also performed a fit using 
the Band function ([Band et al. Ill993l ). The resulting parameters are listed in Table [2] as fit "2". 



Interval "b" is best represented by a Band function with the parameters displayed in Table [2] 
as fit "3" . To estimate the significance of the spectral evolution from interval "a" to interval "b" we 
fit the count spectrum in interval "a" assuming a Band function with a fixed high-energy spectral 
index obtained by the best fit of interval "b" (/3 = -2.1), and we estimate the number of expected 
events in the LAT detector. Based on several realizations, the average number of expected events 
in interval "a" is 2.6, for a chance probability of observing zero counts of about 7%. Therefore, 
we can only suggest that a spectral evolution characterizes the temporal behavior of this burst, 
somewhat similar to that ob served in the first portion of the emission of the long, bright GRB 
080916C (|Abdo et al.ll2nn9cl v Fit "4" shows the result of the fit when a COMPT model plus a 
power-law is adopted. Even though the fit is not favored, and the statistics are very low for any 
conclusive remark, this model has some interesting implications, such as the possible presence of 
an extra component in LAT data. 

The last entry in the table is related to interval "c" which is best represented by a simple 
power-law. Fig. [2] shows the count spectrum for the Fit "3" and Fig. [3] the stacked i^Fi, plot for Fit 
1,3,4 and 5, where the 68% CL are computed from the covariance matrix provided by the fitting 
routine (rmfit). 

Systematic errors are due to uncertainties in the effective areas of the different detectors, 
energy resolution and background estimation. The most important contribution arises from the 
uncertaintie s related to the effe ctive areas. For the LAT they have been derived from a study of the 
Vela pulsar dAbdo et al.ll2009al ) and are 10% below 100 MeV, 5% around 1 GeV and 20% above 10 
GeV. We adopt a 10% uncertainty in the Nal and BGO effective area (both overall normalization 
and slope). We propagate the systematic errors on the effective areas to the parameters of the 
model, and we compute the systematic errors associated to the flux/fluence values. Systematic 
errors in each case are comparable or smaller than the statistical errors quoted in Table [2j 



5. Discussion and interpretation 



GRB 081 024B was the first short GRB with observe d emission above GeV energies. Unlike 



GRB 930131 (Kouveliotou et al. 



1994 



Sommer et al.lll994l '). a 3 GeV photon from GRB 081024B is 



well correlated with the second low-energy pulse. The cross-correlation function of the light curves 
between 30-100 keV and 100-300 keV shows no strong signature of spectr al lag larger than 30 m s, 
which is consistent with the negligible spectral lags in other short GRBs (iNorris &: Bonnellll2006l ). 



While the majority of long GRB spectra are well fitted by the conventional Band func- 
tion, previous spectra l analyses of short GRBs have mostly used the cutoff power-law function 
(jGhirlanda et al.ll2004l : iMazets et al.ll2004l ). The exponential cut-off implies that the bulk motion 
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of short GRBs is not necess arily ultra-relat ivistic, owing to the compactness problem for high- 
energy photons above m-eC^ ( Meszarod |2002| ) . which becomes less severe (see ( Nakarl 120071 ) for a 
quantitative estimate of the Fmin in this case). This difference between long and short GRBs may 
be due to poor counting statistics at high energies in short GRBs, stressing the need for a larger 
sample with sufficient high-energy photons in MeV-GeV bands. 

The delayed onset of a GeV pulse, which is frequently found in other LAT-detected bursts such 
as GR B 080916C or GRB 080825C, may be explained by the different physical pa rameters for two 
pulses (jAbdo et al.ll2009d ). 77 pair-production opacity eff ect (iGranot et al.l 120081). or acceleration 



timescale of high-energy protons for hadronic r nodels (e.g 



Razzaque et al.ll2005l : iDermer fc AtovanI l2006l : lAsano Sz Inoue 



Rachen &: Meszaroi 
2OO7I : 



1998 



Asano et al 



Dermer 2002 



20091 ). The long- 



lasting tail of GeV emission is also common to GRB 080916C, GRB 080825C and GRB 090510. 
One possible interpretation is that the long tail is synch rotron or synchrotron self-Compton (SSC) 
emission during the afterglow phase (iGhirlanda et al. Alternatively, the GeV afterglow 

emissi on may originate from cas cades induced by ultrarelativistic hadrons accelerated by the blast 
wave (iBottcher &: Dermeiill998l ). In these afterglow scenarios even the 3 GeV photon at t ~ 0.5s 
could have an af terglow origin, and t he delayed onset of GeV emission is also naturally explained. 
The onset time (jMolinari et al.l 120071 ) and the hard spectrum for interval c do not contradict the 
afterglow scenarios. Early afterglow models for this l ong lasting tail with synchrotron emission 
(iHe &: Wan g 200!^) and SSC emission (ICorsi et al.ll2009l ) are actually proposed. On the other hand. 



Corsi et al.l (|2009l ) pointed out that SSC emission from late internal shocks can be an alternative 



interpretation for the long lasting tail. 

Finall y, in Fig. [H we show the high-energy fluence versus low-energy fluence diagram, as sug- 
gested by (jLe fc Dermeii l2009l ) . comparing 6 LAT GRBs. This is considered as a way to search 
for separate classes of GRBs and, specifically, spectral differences between the short-hard and long 
duration GRB classes. The sample is limited and therefore no definite conclusion can be drawn, 
but we note that the two short bursts so far detected by the LAT occupy a region where the en- 
ergy emitted at high energy is greater than the energy emitted at lower energy (with a ratio ^ 1) 
suggesting that short GRBs may have an higher efficiency in emitting 7-rays. Nowadays only two 
short GRBs have been detected with a significant emission at GeV energies; we expect that in the 
next decade we will be able to significantly increase the Fermi short GRB sample and understand 
their similarities and differences from the long, soft GRBs. 



The Fermi LAT Collaboration acknowledges support from a number of agencies and institutes 
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Research Council and the National Space Board in Sweden. Additional support from INAF in Italy 
and CNES in France for science analysis during the operations phase is also gratefully acknowledged. 
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Fig. 1. — Multi-instrument light curve for GRB 081024B. The top panel shows the sum of the 
background subtracted signal from two Nal detectors. The second panel is one BGO detector. 
The third panel shows all the events recorded by the LAT, without any selection on the quality of 
the events (background subtracted). The fourth panel shows the selected "transient" events above 
100 MeV. The energy of events is reported at the right axis of the plot. The arrival times, the 
reconstructed positions and the energies of the 11 events are reported in Table [H 
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Fig. 2. — Count spectrum in interval "b", fit "3" . The data points are the rates from GBM (red 
circles: Nalg, green squares: Nalg, yellow crosses: BGOi) and LAT (blue diamonds). The predicted 
rates in the various detectors are obtained by folding the best fit model (Band function) with the 
response of the detectors, and are displayed as continuous cyan lines. 
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Fig. 3. — Analytical representation of the source photon spectra. From top to bottom: fit "1" in 
interval "a" (COMPT model); fit "3" in interval "b" (Band); fit "4" in interval "b" (COMPT plus 
power-law); fit "5" in interval "c" (power- law). 
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Fig. 4. — Fluence-fluence diagram for 6 LAT GRBs. Black squares are classified as short GRBs 
and filled grey circles represent long GRBs. The dashed lines indicate the 1, 0.1, 0.01 values of the 
ratios between the 100 MeV- 10 GeV fluence and the 20 keV-2 MeV fluence. 
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Table 1. Selected LAT "transient" events with energy above 100 MeV. 





Time-Ti^ig 
s 


Energy'^ 


RA 


Dec 


PSF^ 


Dist.4 




MeV 


degrees 


degrees 


degrees 


in PSF 


1 


0.229 


145 ± 22 


325.3 


21.4 


2.6 


0.8 


2 


0.248 


101 ± 17 


323.2 


18.7 


3.4 


0.7 


3 


0.320 


442 ± 44 


322.6 


20.6 


1.1 


0.6 


4 


0.342 


140 ± 21 


322.2 


20.2 


2.7 


0.4 


5 


0.391 


441 ± 44 


323.6 


20.9 


1.1 


0.7 


6 


0.406 


368 ± 41 


323.1 


18.8 


1.2 


1.2 


7 


0.551 


3070 ± 230 


322.9 


21.2 


0.2 


0.2 


8 


1.223 


350 ± 39 


324.1 


20.7 


1.3 


1.0 


9 


1.986 


143 ± 22 


325.0 


24.2 


2.6 


1.4 


10 


2.184 


1680 ± 130 


322.7 


21.5 


0.4 


1.0 


11 


2.801 


386 ± 41 


322.6 


20.1 


1.2 


1.0 



Note. — ^ Arrival time with respect to the GBM trigger time Tt^ig = 
246576161.864; ^ reconstructed energy with estimated error. Errors are esti- 
mated using Montacarlo simulations from the width of the reconstructed energy 
distribution; evaluated 68% containment radius from the PSF for "transient" 
events.'* distance from the GRB position, in PSF units. 



Table 2. Fit parameters for the time-resolved spectral fits. 



Fit Interval Detectors Model Epeak a /3 C-Stat/DOF Fluence Fluence 

(20 kcV - 2 MeV) (100 MeV - 10 GeV) 







(s) 






(MeV) 








erg/cm-^ 




erg/cm^ 


1 

2 


"a" 
"a" 


, 0-0.1 
, 0-0.1 


Nal+BGO 
Nal+BGO 


COMPT 
BAND 


9 7+4.3 
^- ' -1.5 

2.8 


-1 03 +"-23 
^•"■^ -0.19 

-1 03+"'23 


< -1.7t 


330/352 
330/351 


(1.7 + 0.3) X 10" 
(1.7 + 0.4) X 10" 


-7 
-7 


< 4.7 X 10^1" 

< 7.3 X 10-8 


3 
4 


"b" 
"b" 


, 0.1-0.8 
, 0.1-0.8 


Nal+BGO+LAT 
Nal+BGO+LAT 


BAND 
COMPT + 
POW 


2Q+1.U 

-1.0 


-1 03 +"■-'''' 
^•'^■^ -0.14 

-0 7+° * 


-2 10 +"-^1 
^■-'-'J -0.14 

-1.OO-0.06 


383/359 

384/358 


(7.9+1.9) X 10" 
(8.8 + 2.0) X 10" 


-7 

-7 


(8.9 + 3.0) X 10"'^ 
(1.1 + 0.4) X 10"*^ 


5 


"c" 


, 0.8-2.9 


Nal+BGO+LAT 


POW 






-1 6+"'' 


274/361 


(4.3 + 3.2) X 10" 


-8 


(4.0 + 2.4) X lO"'' 



Note. — For each time bin ("a", "b" and "c"), and for ea<;h fit (1 to 5), we report the detectors used in the fit, the name of the model 
used, and the value of each spectral parameter. Asymmetric errors are obtained from the profile of the Cash statistic. The pivot energy of the 
power-law function is set to 100 MeV. t The upper limit for /3 has been obtained looking at the profile of the Cash statistic. Although in this fit 
/3 is basically undetermined, we report it for reference. 



